Topological Insulators with Commensurate Antiferromagnetism 
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We study the topological features of non-interacting insulators subject to an antiferromangetic 
(AFM) Zeeman field, or AFM insulators, the period of which is commensurate with the lattice 
period. These insulators can be classified by the presence/absence of an emergent anti-unitary 
symmetry: the combined operation of time-reversal and a lattice translation by vector D. For 
AFM insulators that preserve this combined symmetry, regardless of any details in lattice structure 
or magnetic structure, we show that (i) there is a new type of Kramers' degeneracy protected by 
the combined symmetry; (ii) a new Z2 index may be defined for 3D AFM insulators, but not for 
those in lower dimensions and (iii) in 3D AFM insulators with a non-trivial Z2 index, there are odd 
number of gapless surface modes if and only if the surface termination also preserves the combined 
symmetry, but the dispersion of surface states becomes highly anisotropic if the AFM propagation 
vector becomes small compared with the reciprocal lattice vectors. We numerically demonstrate the 
theory by calculating the spectral weight of the surface states of a 3D TI in the presence of AFM 
fields with different propagation vectors, which may be observed by ARPES in Bi2Se3 or FfeTes 
with induced antiferromagnetism. 



I. INTRODUCTION 

The inception point in the study of topological phases 
in condensed matter systems is the integer quantum Hall 
effect (IQHE). Since the Hall conductivity is odd under 
time-reversal symmetry (TRS), the topologically non- 
trivial states will only occur when TRS is broken. The 
connection between the observed quantized Hall conduc- 
tance and topology is made through the calculation of 
the Chern number which, in the context of the IQHE, 
corresponds to the number of chiral edge states present 
and is a global quantity contributed by all of the oc- 
cupied states in the system 1 . Nevertheless, the vast 
majority of topological band insulator states that have 
been experimentally confirmed 2-5 , are stabilized by the 
presence of TRS in which strong spin-orbit interactions 
invert the bands which allows the creation of a mass 
boundary 6-8 that traps surface or edge states at the in- 
terface of two bulk band insulators or that of a bulk in- 
sulator and vacuum, commonly referred to as topolog- 
ical insulators (TI). While the theory associated with 
TRS TI is well understood, the presence of TRS in a 
material significantly constrains the number of potential 
topological materials. Therefore, it is important to go 
beyond the well-established paradigm of non-interacting 
time-reversal symmetric TI to look for novel topologi- 
cal phases. To do so, one must generalize the theory of 
topological classification to encompass a wider spectrum 
of materials 9-1 . 

Of the materials not covered in the standard theory 
of TRS TI, one of the most interesting classes are mag- 
netic materials. The established theory for TRS insula- 
tors cannot be applied to classify any magnetic materials 
as their ground states break TRS. Ferromagnetically or- 



dered materials have been explored theoretically as their 
broken symmetry leads to the the possible condensed 
matter realization of Wcyl fcrmions 16 ' 20 ' 21 . (While the 
magnetic structure studied in Rcf . [20] has zero total mag- 
netization, we still refer to it as fcrromagnctism because 
of the unbroken translational symmetry.) Yet, to this 
point, there has been comparatively little work on anti- 
ferromagnetically ordered materials. In a work by Mong, 
Essin and Moore 22 , the insulators subject to an AFM 
Zeeman field of propagation vectors (0, 0, it), (tt. tt, 0) and 
(tt, 7t, tt) are studied. The spin structure in these insula- 
tors can be seen as sets of ferromagnetically ordered 2D 
planes of spins with alternating magnetization, stacked 
along [001], [011] and [111], respectively. Mong et al show 
that a new Z2-indcx can be defined for 3D AFM insula- 
tors of this type. Nevertheless, real magnetic materials 
usually exhibit more complex magnetic structures, e.g., 
having non-collinear configurations as spiral magnetism, 
having a large magnetic unit cell (small propagation vec- 
tors) or having more than one propagation vector. These 
possibilities are specially relevant for AFM field induced 
by dopants in an otherwise non-magnetic insulator like 
Bi 2 Se 3 or Bi 2 Tc3. Therefore, it is necessary to build a Z 2 
classification for AFM insulators with generic magnetic 
ordering, which we do in this paper. 

Let us recall that the existence of Kramers' 
degeneracy 23 is necessary for defining the Z2 index for 
TRS insulators 24-28 . Kramers' degeneracy is the exact 
double degeneracy of every single particle state in a spin- 
1/2 system with TRS. We begin by studying the breaking 
or preservation of Kramers' degeneracy in AFM insula- 
tors. While it is clear that ferromagnetism always breaks 
a Kramers' pair into energetically separated spin- up and 
spin-down states, we show that AFM preserves the de- 
generacy of almost all single particle states, with excep- 



tions detailed in later sections, if there is a lattice trans- 
lation vector D that inverts the spin at every site. This is 
because if such D exists, we can define a new anti-unitary 
operator 22 @s = @Tb, where is the time-reversal sym- 
metry operator, and Td represents a translation by D. 
@s is a symmetry because both and Td invert the 
magnetization and so their combination recovers it. It is 
this symmetry that ensures that Kramers' degeneracies 
are preserved for almost all single particle states. AFM 
insulators that are invariant under ©s will be referred to 
as 05-symmetric insulators hereafter. 

We may arrive at the same result beginning from a 
different perspective, by treating the AFM field, a static 
field coherently coupled to the electrons in the form of 
c£(r)(J(r) • cr) TT 'C T '(r), where a is the spin operator and 
t, t' the spin indices, as perturbation and expressing the 
single-particle Green's function to an arbitrary order of 
perturbation. This allows us to calculate the poles of the 
Green's function to obtain the energy spectrum. From 
this perturbative point of view, the energy spectrum of 
the TRS insulator is modified by the self-energy contri- 
bution from scattering by the AFM field, and Kramers' 
degeneracy would require that the self-energy term be 
invariant under TRS. The full self-energy counts in all 
scattering processes that send the particle to its initial 
state. Since the AFM field changes sign under TRS, 
the scattering contribution from even number of scat- 
terings are invariant under TRS and those from odd 
number of scatterings are variant. Therefore, the total 
self-energy term is invariant under TRS if and only if 
there are only even number of scatterings. Any com- 
mensurate AFM field decomposes into a finite series of 
normal modes, namely J(r) = J2i=i m M(Qi)e lQi ' r , 
where Qi=i ) ..., m are the set of propagation vectors of 
the AFM field. Consider an electron with initial mo- 
mentum k, then after N scatterings its momentum be- 
comes k' = k + Z1Q1 + Z2Q2 + •■■ + z m Q m , where 
z\ + z-2 + ... + z m = N. If this scattering process con- 
tributes to the self-energy, we have k' = k + G, where 
G is any reciprocal lattice vector. Preservation of the 
Kramers' degeneracy thus requires that only N = even 
terms are nonzero. From this, we obtain the necessary 
and sufficient condition for existence of a Kramers' de- 
generacy: for any set of m integers (z\, ..., z m ) satisfying 
Y^i z iQi = G, there is J2 i Zi = even. Mathematically, 
satisfying this condition is completely equivalent to the 
existence of a lattice translation vector D that inverts the 
AFM field. This explicitly indicates that the propagation 
vectors alone determine whether a Kramers' degeneracy 
of 0s is preserved or broken, while any other detail in the 
magnetic structure is irrelevant. Also, we use the pertur- 
bation theory to estimate the split of the Kramers' pair 
when the above condition is violated. 

Having established the necessary and sufficient con- 
ditions for the Kramers' degeneracy in a commensurate 
AFM, we further investigate if this allows us to define a 
Zi topological invariant. It is known that one can define 
the Zi invariant using the Pfaffians of the sewing ma- 



trix of time-reversal operator at time-reversal invariant 
momenta (TRIM) in 2D and 3D TRS insulators 24 ' 26 ' 27 . 
We prove that in ©s-symmetric AFM insulators, a new 
^-invariant can be defined only in 3D (or certain higher 
dimensions) using Pfaffians of the sewing matrix associ- 
ated with S at half of the TRIM, while in 2D and ID, 
the same definition gives a gauge variant quantity. This 
is because the d-dimcnsional AFM Z2 number is defined 
in its d — 1-dimcnsional subsystem that belongs to the 
symplectic class (All in the A-Z classification 29 of all 
Hamiltonians) , and OD and ID symplectic Hamiltonians 
are always topologically trivial 9-11 . However, when in 
addition to ©s there is also spatial inversion symmetry, 
one can use inversion eigenvalues at half of the occupied 
bands to define a new ^-invariant in arbitrary dimen- 
sions. 

In TRS insulators there is a correspondence between 
the bulk .^-invariant and the existence of gapless bound- 
ary modes 30 : when the Z2 number is non-trivial, there 
are odd number of gapless (Dirac) modes at the boundary 
of the system. Here we show that the existence of gap- 
less surface states in a 3D AFM insulator requires two 
conditions: (i) it has a non-trivial ^-invariant in the 
bulk and (ii) the surface termination also preserves the 
©s-symmetry, i.e., the ©s symmetry is unbroken on the 
surface. If the AFM is small, one may say that the Dirac 
surface modes of 3D TRS TI arc preserved in presence of 
a ©s-symmetric AFM. However, in a 2D ©s-symmetric 
AFM insulator, the surface modes are not protected from 
being pushed into the bulk by tuning the surface chem- 
ical potential. This is in exact agreement with the fact 
that the new ^-invariant is only well-defined in 3D Q$- 
symmetric AFM insulators and not in lower dimensions. 

While a TI that has intrinsic AFM ordering is yet to 
be experimentally established, the theory can be applied 
to TI thin films with induced AFM ordering by an AFM 
substrate. In principle one can use the substrates to in- 
duce various AFM orderings in the TI thin film, but lat- 
tice matching at the interface is a practical issue. The 
lattice constants of the substrate should be nearly identi- 
cal to those of the TI at the interface, or at least commen- 
surate with them. If the AFM substrate is ©s-symmetric 
in the redefined basis, the Dirac point on the surface of 
the TI would be unbroken by the induced AFM, while the 
surface state dispersion becomes anisotropic along the di- 
rections parallel and perpendicular the AFM propagation 
vector. 

It is straightforward to extend the discussion to topo- 
logical superconductors (TSC) with coexisting AFM or- 
ders. Given a d-dimensional TSC in the Dili class (SC's 
with non- negligible spin orbital coupling and with TRS), 
which has an induced or intrinsic AFM order, if the AFM 
ordering is ©s-symmetric, then we show that the AFM 
TSC has the same classification as a d — 1-dimcnsional 
TSC without AFM. Therefore AFM SC's in 3D and 2D 
have Z2 topological classifications 9 ' 11,31 . 

The paper is organized as follows. In Sec. II, we show 
that a new type of degeneracy is preserved in an AFM 



insulator if and only if it is ©s-symmetric. In Sec. Ill, 
we use perturbation expansion of self-energy in the sin- 
gle particle Green's function to re-derive the condition for 
Kramers' degeneracy, and, when the condition is not met, 
the relation between the energy splitting of a Kramers' 
degeneracy by the AFM and the propagation vectors of 
the AFM. In Sec. IV, the results obtained in the previ- 
ous sections are used to classify 3D AFM insulators by 
a ^-invariant. In Sec.V we study the boundary modes 
of 2D and 3D AFM insulators, reaching the conclusion 
that only 3D Os-symmetric AFM insulators may have 
protected gapless Dirac modes, yet anisotropic, on the 
surface that is also Os-symmetric. In Sec. VI, we discuss 
how the theory may be applied in real materials and also 
how an analysis similar to what we perform in AFM insu- 
lators can be easily extended to AFM superconductors. 
We conclude the work in Sec. VII. 



II. KRAMERS' DEGENERACY IN AN AFM 
INSULATOR: PROOF BY SYMMETRIES 

Throughout the paper, we assume that the electrons 
in the AFM insulator can be modeled by a TRS tight- 
binding Hamiltonian that coherently couple to an anti- 
ferromagnetic Zeeman field: 



H = H Q + H M , (1) 



where 



H = J2 K T ' PT '(k)cUk)cp T ,(k) (2) 

k,a,/3,r,r' 

is the hopping part that is TRS and 

H M = J2 (J(rWrrOcUr)w(r) (3) 



is the part coupled to the Zeeman field J(r). In the 
above expressions and hereafter, Greek letters denote the 
orbital degrees of freedom inside a unit cell with excep- 
tions of r, t' which refer to up and down spin compo- 
nents. We assume ai=i <j's to be the basis vectors of 
the d-dimensional lattice in which the model is embed- 
ded; and bi=i v .. )( j's are the corresponding reciprocal lat- 
tice basis vectors. The AFM field J(r) breaks the lattice 
translation symmetry, but as long as its periods are com- 
mensurate with lattice vectors, it enlarges the unit cell 
of the crystal lattice. The enlarged unit cell, or magnetic 
unit cell, is given by the basis vectors a^ s, which is in 
general a linear combination of the a^ 's with integer co- 
efficients. The counterpart of the enlarged unit cell in 
the reciprocal space is the folded Brillouin zone (BZ), or 
the magnetic BZ (MBZ) spanned by basis vectors bV's 



2irSi 



satisfying af 1 ■ h 1 } 1 — *,,u tJ . 

Before proceeding to discuss the degeneracy in a time- 
reversal breaking AFM system, it is helpful to briefly 
review the concept of Kramers' degeneracy in a TRS sys- 
tem, and the key role it played in the ^-classification of 



TRS insulators in 2D and 3D. The time-reversal operator 
in the single-fermion sector of Hilbcrt space, denoted by 
is defined as: 



= K{ia y 



(4) 



where K is the complex conjugate operator. From this 
definition one can see that is antiunitary and squares to 
minus the identity. These two properties guarantee that 
each single particle state is at least doubly degenerate 23 . 
This degeneracy protected by time-reversal symmetry is 
called Kramers' degeneracy and the degenerate doublet 
referred to as a Kramers' pair. In a translationally in- 
variant system, one can further show that an eigenstatc 
with momentum k is sent by to an equal energy eigen- 
state at momentum — k. Specially, at center or corner of 
BZ where we have — kmv = kj n „ mod G hence every en- 
ergy level of H(k inv ) must be doubly degenerate. These 
points at BZ center and corners are called time-reversal 
invariant momenta (TRIM) and the Kramers' degeneracy 
at TRIM make possible the definition of a ^-invariant 
in 2D and 3D TRS insulators. In order to formally de- 
fine this invariant, one uses the sewing matrix, which is 
defined as 



B mn (k) = (Va-k)|0|Vv>(k)), 



(5) 



where |^> m (k)) is the energy eigenstate at k in the m-th 
band. The sewing matrix is antisymmetric at all TRIM 
and this antisymmetry enables the definition of the Pfaf- 
fians at TRIM which form the basis of the definition of 
the i^-invariant Sq, given by (Fu-Kane formula 27 ): 



i-i) 5 ° = n 



Pf[g(kin„)] 

y/det[B(k inv )] 



(6) 



It should be noted that although it appears that <5o only 
depends on the band structure at TRIM, it implicitly 
depends on the band structure in the whole BZ, as a 
smooth gauge is required for this definition. 

With the role of Kramers' degeneracy defined in TRS 
insulators, let us turn our attention to the AFM insula- 
tors. Any AFM breaks TRS because TRS reverses all 
spins leaving the orbital and spatial components invari- 
ant as: 



Mm®- 1 = -H 



M- 



(7) 



Therefore, Kramers' degeneracy in its original meaning 
does not exist in any AFM insulator. However, this 
does not necessarily mean that the spectrum is non- 
degenerate. In fact, each energy level may still be doubly 
degenerate, but the two states are not each other's time- 
reversed counterparts. In that case, there must be some 
other symmetry that protects the pair from splitting in 
energy. For a generic AFM TI, the only symmetry is the 
magnetic translation symmetry by any magnetic super- 
lattice vector L A/ (a superposition of af's with integer 
coefficients), i.e., 



[H,Tj 



L"J 



0, 



(8) 



where Tl is the operator for a translation by a lattice 
vector L. Eq.(8) states that one can find a set of basis 
vectors that are common eigenstates of H and T-^m . The 
definition of T^m dictates that its eigenvalues must be 
of the form exp(ik • L M ), where k is constrained within 
the first MBZ. This symmetry gives us a band structure 
defined in MBZ, but does not protect any degeneracy 
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FIG. 1. AFM configurations that are (a,b) ©s-invariant 
and (c) ©s-variant. In (a), we plot the typical spin struc- 
ture of a CuC>2 layer in any parent compound of cuprates 
superconductors 32 . In (b), we show the spin structure of va- 
cancy doped iron-based superconductor Rb2Fe4Ses . In (c), 
we plot the ground state of a classical Heisenberg model on a 
triangular lattice. The basis vectors of the lattice, ai,2, and 
those of the magnetic superlattice, a*^ are also shown. Note 
that in (a) and (b), a translation by D = ai inverts all spins, 
while in (c) one cannot find such a vector. 



In a large class of magnetic structures, examples of 
which arc shown in Fig.l(a,b), there exists a special lat- 
tice vector D, where, after a translation by D, all spins 
flip their signs. Mathematically, we may define this as 



{f u ,H M } = 0. 



(9) 



The translation vector, D, is unique only up to some 
magnetic translation vector L M . We now suppose that 
both Di and D2 satisfy Eq.(9), then we have 



[Td 1 -d 2 ,H m ] — [T Dl T D2 , H M ] — 0, 



(10) 



which indicates Di — D2 = L . Using the translation^ 
invariancc and TRS of the free part Hq 7 we have 

[H Q ,f- D ] = [H o ,e] = 0; (11) 

using Eq.(7) in conjunction with Eq.(9) we have 

[H M ,fnQ] = 0. (12) 

Combining Eq.(ll) and Eq.(12), we obtain a new sym- 
metry of the full Hamiltonian, s = 9Tb: 



[H,Q s } = 0. 



(13) 



Since Tb is unitary and antiunitary, Q$ is antiunitary, 
namely: 



On the other hand square of 0s is 

el = e 2 T 2D = -f, 



2D, 



(14) 



(15) 



where we have used the fact that O and T T commute. 

We now turn our attention to the degeneracy in the 
spectrum of H. Suppose |k,n) is an eigenstate of both 
H and T l m with k a wave vector in the MBZ, then 

H\k,n) = E n (k)\k,n) (16) 

f L M|k,n) =exp(ik-L M )|k,n), (17) 

then from Eq.(13) we have 

HQ s \k,n)=E n (k)Q s \k,n), (18) 

meaning that 65 |k, n) is also an eigenstate of H with the 
same eigenvalue. To check if |k, n) and Os|k, n) are same 
or different state, we calculate the overlap using Eq.(14) 
and Eq.(15) as 



(k,n\Q s \k,n) = (G 2 s k,n\Q s \k, 



(19) 



(k,n|T_ 2D 9 s |k,n) 



At this point we notice that 2D is a magnetic superlattice 
vector as it leaves Hm invariant, which implies that 



-2ikD 



(k, n|T_ 2 D = (k, n\e 
Substitute Eq.(20) into Eq.(19) and have 



(k,n\e s \k,n) 



-2ik D 



k,n|8g|k,n). 



(20) 



(21) 



Eq.(2I) shows that unless 2k • D — (2n + l)ir, we have 
(k, n|0s|k, n) = 0, ensuring that |k, n) and 05|k,n) are 
two orthogonal states. Beyond this, it is crucial to un- 
derstand how the state Os|k, n) transforms under a su- 
perlattice translation or 



T-lm & s \k,n) 



= 6 s T L Af|k,n) 
|k,n) 



(22) 



9 s e lkL 

„-ikL A ' 



9 s |k,n) 



Therefore 65 |k, n) must be proportional to a state with 
wavevector — k, or mathematically 

|-k,n)= Y, B mn {k)\k,n), (23) 

n£occ 

where m,n are band indices (to 7^ n in general due to 
degeneracy) and B mn (k) is called the sewing matrix of 
O5 at k. Moreover, from Eq.(22), we also notice that 
other than at TRIM 8s|k, n) is orthogonal to |k, n) for 
a generic k because they have different eigenvalues under 
superlattice translation. 

To this point, we have shown that at all non-TRIM, the 
eigenstate at k must be degenerate with another state at 
— k. Furthermore, at TRIM for which exp(z2ki nl ,-D) = 1, 
each level at ki nv must be doubly degenerate and that 
at TRIM for which exp(i2ki„„ • D) = —I, the levels at 
ki nv are gcncrically non- degenerate if there is no other 
symmetry present in the system. As we know that the 
TRIM with degenerate energy levels are important for 



the definition of ^-invariant in a TRS insuiator, thus it 
is important to distinguish the TRIM with exp(iZk inv ■ 
D) = — 1. Since 2D is a supcrlattice vector, we have the 
general form of D as 



D= £ x ia f/2 + L M , 



(24) 



=l,...,d 



where xi is either zero or one while there is at least one 
xi = 1. On the other hand, all TRIM can be represented 
by 



k in „= Y V&i*/ 2 ' 



(25) 



where yt is either zero or one. Then if exp(i2\ii nv • D) = 
— 1, we have, using af f • bj' 1 = 2n5ij, 



Y. XiVi e odd. 



(26) 



-l,...,d 



There are 2 d TRIM, and for a given set of {x{\, exactly 
half of them satisfy Eq.(26). The proof goes as follows. 
Since there must be one Xi that is nonzero, we can as- 
sume x\ — 1 without loss of generality. Then for any set 
of {y 1 ,y 2 ,---, yd} that satisfies Eq. (26), {1 -3/1,2/2, -,2/d} 
satisfies J^ Xiyi = even. So there is a one-to-one cor- 
respondence between TRIM that satisfy Eq.(26) (called 
from now on A-TRIM) and those that do not (called 
B-TRIM). At an A-TRIM, each level is non-degenerate 
while at a B-TRIM, each level is doubly-degenerate. 
For a simple example, consider a 2D antiferromagnet 
with a single propagation vector q = (x,ir) (Fig. 1(a)). 
This means bf = (tt/2, tt/2) and hf = (tt/2, -tt/2), 
and therefore A-TRIM are (tt/4, ±tt/4) and B-TRIM are 
(0,0), which is always a B-TRIM, and (tt/2,0) corre- 
sponding to (2/1,2/2) = (1, !)■ 



III. KRAMERS' DEGENERACY IN AN AFM 
INSULATOR: PROOF BY PERTURBATION 

In the previous section, we have proved, by using only 
properties of symmetry operators, that a new type of 
Kramers' degeneracy is protected by a combined symme- 
try 0s . The cases of single propagation vectors such as 
(0,0, 7r), (71-, 7T, 0) or (7r,7r,7r) have been shown in Rcf.22, 
but here we extended it to arbitrary commensurate AFM. 
The validity of the statement is independent of the partic- 
ular strength of AFM or the specific magnetic structure. 
We here give a more intuitive approach to understand- 
ing this result. We start from another perspective by 

I 



treating the AFM as a perturbation field applied to an 
otherwise TRS system, and study how the perturbation 
breaks/preserves the Kramers' degeneracy of the TRS 
system. 

The single particle Hamiltonian of Eq.(l) can be 

-q+G qi q 2 -qi - q2+G 



(a) q 

£(<u,k) = 



k + q 



k + qi k + q t + q 2 

-qi - qa - q 3 +G 



q3 



(b) 



k + qi k + qi + q 2 k + qi + q 2 + q 3 
(c) 



1 



uj - H a (p) 



= M(p) a 



FIG. 2. (a) is the Feynman diagram expansion of the self- 
energy induced by the scattering of AFM field, (b) represents 
the free electron propagator in the momentum space without 
the AFM field, (c) is the vertex representing the scattering 
by the AFM field of momentum p. 



rewritten as 

H= Y, ?W.'(k)cUk) C/3T ,(k) (27) 

k,a,/3,r,r' 

+ Y (M(Q,)-cr TT /)4 r (k+Q i )c^(k). 



where M(Qj)'s are the Fourier components of the Zee- 
man field J(r). The energy levels at k present themselves 
as the poles of the single-particle Green's function: 



G(w,k) 



1 



w-'Ho(k)-E(w,k)' 



(28) 



where E(w, k) is the self-energy. When the AFM field 
is weak, the self-energy can be expanded in terms of the 
number of scatterings of an electron whose initial and 
final momentum are k plus a reciprocal lattice vector 
G. This scattering process can be visualized using Feyn- 
man diagrams, as shown in Fig. 2(a), where the first three 
terms represent the electron being scattered by the AFM 
field once, twice and three times before returning to the 
initial state. Applying the conventional Feynman rules 
given in Fig.2(b,c), the self-energy is given by perturba- 
tion series: 



E(w,k) = Ei(o;,k) + S a (w,k)- 
E„(w,k)= Y 

q I e{Qi,...,Q m },qi + ..+q„=G 



+ E n (w,k) + ... 

[M(qi)-o-]G (u;,k + qi) 



[M(q„ 



<7]G (u;,k + qi 



(29) 



•q„_i)[M(q n ) -a] 



where the free Green's function 



Go(w,k) 



-flo(k)' 



(30) 



We note that for the Green's function, w should have 
an infinitesimal imaginary part, but since we are only in- 

I 



tcrestcd in the position of real poles, we take w to be a 
real number. Due to the TRS of the unperturbed Hamil- 
tonian, all poles of the free Green's function Go(us, k) arc 
doubly degenerate. If there is a splitting of the degen- 
eracy, it must be due to the self-energy term E(w,k). 
Though in general S(w, k) is not hermitian, it is hermi- 
tian in this problem as only elastic scattering is present 
(and using M(Q 4 ) = M*(-Q 4 )): 



Et( w ,k)= Y, [M(-q„)-a]G (a;,kH-q 1 + ... + q„_i)...[M(-q 2 )-a]Go(a;,k + q 1 )[M(-q 1 )- C 7] 

qi G{Qi.....Q m },qi + ..+q„=G 

Y [M(-q„) • a]G {oJ, k - q„)...[M(-q 2 ) • a]G (co,k- q„ - ... - q 2 )[M(-qi) • a] 

- q< e{Qi,...,Q m },-q„-...-qi=G 

= E n (w,k). (31) 



Therefore all poles of G(ui, k), like those of Gq(uj, k), are 
still real. 

We will now focus on the symmetry property of each 
term in the expansion of E(w,k). Under time-reversal 
operator of Eq.(4), the free Green's function transforms 
as 

eGo(^,k)e- 1 =G ( W ,-k). (32) 

The AFM vertex from the expansion transforms as 

9M(Qi) • aQ- 1 = -M(-Qi) ■ a. (33) 

Therefore the n-th order self-energy transforms according 
to 



es n ( w ,k)e- 1 = (-i) ri s n ( w ,-k). 



(34) 



Since for n £ even, 0E„(w,k)0 l = £ ra (u>, — k), then if 
S„ = for any n € odd, we have 



9E(w,k)9- 1 = E(w,-k). 



(35) 



Using Eq.(35) in conjunction with Eq.(4), we know for 
each eigenstate |u(w,k)} of "Ho(k) + E(w,k), 0|M(w,k)) 
must be an orthogonal eigenstate of Ho(— k) + E(w, — k) 
with the same eigenvalue. Therefore, as far asn£ odd 
terms are vanishing in the self-energy expansion, the 
AFM as a perturbation will not split the Kramers' de- 
generacy in a TRS system. Physically, this can be easily 
understood as follows: if an electron must experience an 
even number of scatterings by the AFM field before go- 
ing back to its initial state, the scattering amplitude does 
not change under time-reversal which changes the sign 
of AFM; therefore the single electron propagator cannot 
'see' the breaking of TRS and its poles remain doubly 
degenerate as in the case where TRS is preserved. 

We now determine the condition under which the odd 
terms in the self energy disappear. Let us start from an 
AFM of a simple type in a ID insulator to gain some 
intuition. We assume that there is only one wavevector 



I 

p = 2n/s, where s > 1 is an integer. In the n-th (n € 
odd) order term of E(w,k), each g, must take either p 
or — p, but because of the constraint J^ q{ — G, we can 
easily see that if s = even, it is impossible to satisfy 
the constraint with an odd number of q^s, or in other 
words, E neo dd(w, k) = 0. Moreover, if s <G odd, there 
are two cases: The first case is when n < s where it is 
still not possible to satisfy the constraint if n € odd, i.e., 
^neodd,n<s(w, k) = 0. The second case is when n > s, 
in this situation one may choose: q± = ... = q s = p and 
Qs+k = {— 1) p for < k < n — s (since n — s = even 
the sum over q s +k vanishes), to satisfy the constraint if 
?i G odd, i.e., E„ Godd! „> s (w,k) ^ 0, thereby breaking 
the degeneracy. Thus, by use of this simple example, we 
know that if s = even, all odd terms in E(w,k) must 
vanish and if s G odd, the lowest nonvanishing odd order 
is the s-th order. 

A similar analysis can be made for the general case. 
In a general commensurate AFM, there arc a finite num- 
ber of propagation vectors denoted by Qi,...,Q m - At 
each scattering, there is q.; € {Qi,...,Q m }, so a generic 
scattering process that contributes to the self energy 
has Zi times of momentum transfer Q^ and satisfies 
5Zj ZiQi = G. Therefore, if we want all odd orders in 
E(o;,k) to vanish, we require that for any set of integers 

Zj=i m satisfying £\ z^Qj = G, J2i z t = even. This is 

the sufficient and necessary condition under which every 
pair of doubly degenerate poles of Gq(io, k) remains de- 
generate in the poles of G(w, k), given that the AFM field 
is small enough to be treated as perturbation. (Unlike in 
the simple ID example, however, when the condition is 
not satisfied, the general formula for the lowest order of 
scattering that breaks the degeneracy is unknown to us.) 

In this and the previous section, we have established 
two equivalent sufficient and necessary conditions for 
the preservation of Kramers' degeneracy in an AFM TI, 
which are: 

(i) there exists a lattice vector D, the translation by 



which flips all spins in the AFM and 

(ii) for any m integers satisfying £\ ZiQi 
J2i z i = even. 



G, then 



The actual proof of their equivalence is purely mathe- 
matical and is hence deferred to the Appendix. Here we 
comment that condition (i) is a direct extension of the 
case discussed in the earlier work by Mong et al 22 , where 
AFM field with propagation vector (0,0, 7r), (0,tt,tt) or 
(7r,7r,7r) satisfies this condition. Condition (ii) is a con- 

I 



straint on the propagation vectors and these vectors 
alone, making explicit the fact that the preservation of 
Kramers' degeneracy has nothing to do with the magni- 
tudes or directions of M(Qi). 

Besides relating the propagation vectors to the degen- 
eracy in the energy spectrum, the analysis of perturba- 
tion helps estimate the magnitude of the energy splitting 
if condition (i) or (ii) is violated. To see this, we re- 
turn to the simplest example of ID AFM with wavevec- 
tor p = 2t:/s. We have mentioned that when s G odd, 
the lowest order of self-energy that breaks the Kramers' 
degeneracy is the s-th order: 



E a («, k) = (M(p) ■ (t)G (u, k+p)(M(p) ■ <t)G (oj, k + 2p)...(M(p) ■ a)G (w, k + (s - l)j>)(M(p) • er) 



(36) 



The magnitude of this term is estimated as: 



(a) 



(b) 



E s (w,fc) 



(37) 



\M\ 



{uj - E{k+p)){uj - E(k + 2p))...(w - E{k + (s - 1») 

\Ml 

~rii=i,...,.-ii^(*)--B(*+*i')r 

To go from the second to the last line, remember that 
since we are interested in the correction to the position 
of the real pole which is at u> — E(k) without AFM, we 
can approximately substitute uj = E(k) into the second 
line. 

From this analysis, we can see that, assuming one can 
treat AFM exchange field perturbativcly, the splitting be- 
comes exponentially small as s, or the period, increases, 
justifying the intuition that when s is large enough, there 
should be no difference between even s (no splitting) and 
odd s (exponentially small splitting). 

One can always use the same analysis to estimate the 
splitting with a more complex magnetic structure. The 
key point is to find the first term in the series expansion 
that has odd number of M(q)'s. Suppose it is the fc-th 
term, then the energy splitting is approximately 



AE ex \M\ 



(38) 



For the simplest example, suppose there are two propa- 
gation vectors satisfying Qi = 2Q 2 , then the lowest odd 
order of scattering that contributes to the self energy is 
three: an electron can gain momentum of Qi in one scat- 
tering and lose momentum Q2 in two scatterings. In this 
case, the energy splitting scales like |M| 3 , independent of 
the actual values of Qi.2- 

To solidify our understanding of the energy splitting, 
we use explicit models of ID insulators with commensu- 
rate AFM. We take a basic four-band model for the TRS 
part of the ID insulator: 

H (k) = (m - cos(fc))r + sin(fc)(r 1 + bT 15 ), (39) 
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FIG. 3. (a)Magnetic structure of a sinusoidal antiferromag- 
net with wavevector Q = 7r/5 described in Eq.(40). (b) The 
energy splitting of the Kramers pair in a ID TRS insulator 
in the presence of AFM shown in (a) at k — as a function 
of the spatial period s of the AFM. Red dots correspond to 
odd and blue empty circles correspond to even s's. (c) Mag- 
netic structure given by Eq.(41) with s = 10. (d) The energy 
splitting of the Kramers pair in a ID TRS insulator in the 
presence of AFM shown in (c) at k — plotted against the 
AFM strength cubed, M 3 . 
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ri = a x ® t 2 , r 2 = 1 

=0,1.2,3 ^» an d Tij = iTiTj. 



9Ty, 13 = 

The time- 



where To 

<y y <8> t x , r 5 = n, 

reversal operator is 6 = AT15. For the AFM part of 
the Hamiltonian, we consider the magnetic structure of 
a sinusoidal spin density wave with the period of s lattice 
spacings (see Fig. 3(a)): 



,, s , ,, 2nx 
J(x) = iWzcos . 



(40) 



Taking the parameters m = 0.5, b = 0.01 and M = 
0.1, we calculate the energy splitting of the Kramers' 
degeneracy at k = in the presence of the AFM field, and 
plot it as a function of n in Fig. 3(b). From this Fig. 3(b), 



it is clear that when s — even, the energy splitting is zero 
(considering the numerical error) and when s £ odd, the 
energy separation decays exponentially with increasing 
s for s ^ 1, as predicted in Eq.(38). Next we consider 
a magnetic structure with two harmonics, as shown in 
Fig.3(c), given by 



must vanish. The Hall conductance can be calculated as 



J(x) = Mz(cos 



2TTX 



1 Attx 
-cos—). 



(41) 



Since Q\ = 2ir/s and Q2 = 47r/s, there is Q\ — 2Q2 = 0. 
This means that an electron at momentum k can return 
to the initial state by three scatterings with the AFM 
held, contributing a TRS odd term to the self-energy 
which is proportional to \M\ 3 . Therefore, according to 
Eq.(38), the energy splitting should be proportional to 
\M\ 3 . In Fig. 3(d), we plot the energy separation against 
\M\ 3 for s = 10 and we find that the result is consistent 
with the our analysis. 

As a final remark, all this perturbative analysis re- 
quires that |M(Qj)| be much smaller compared with the 
energy separation \E„(k) — Eik + Qj)| and hence does 
not apply to flat band Hamiltonians where E n (k) = 
E(k + Qi). The proof by symmetry presented in the 
previous section, however, applies to all cases. 



IV. ^-INVARIANT FOR AFM INSULATORS 

In previous sections, we have derived the conditions 
under which the Kramers' degeneracy can be preserved 
in an AFM insulator. Specially, we have derived that at 
half of the TRIM (called B-TRIM) each energy level is 
doubly degenerate. In a TRS insulator the degeneracy at 
TRIM is key to defining the Zi topological invariant in 
2D and 3D. Therefore, in this section, we seek to define 
the Zi invariant in an AFM insulator while restricting 
the discussion to dimensions d < 3. 

We start by inspecting the sewing matrix associ- 
ated with 85 symmetry at B-TRIM. Using Eq.(14) and 
Eq.(15), we have: 

BmnO^inv) = (4>m 0<-inv) | &S 1 4>n (Knv)) (42) 

= (6!Vn(ki„„)|6,s|V>m(kjmO) 

= -(lpnO<-inv)\&s\lpmO<-inv)) 

This tells us that the sewing matrix is antisymmetric at 
all B-TRIM, which allows us to define the Pfaffian of the 
sewing matrix and, with it, define a Zi quantity <5q : 



(_!)*, 



n 



ki„„eB-TRIM 



Pf[g(kin„)] 

y/det[B{k inv )} ' 



(43) 



For the Zi index in Eq.(43) to be well-defined, we need a 
smooth B(k) to have a consistent sign for the square root. 
This essentially requires that Stoke's theorem applies to 
any non-contractible closed loop in the MBZ, or cquiva- 
lently, that any component of the Hall conductance af? 






2D plane in 3D BZ 



(d l u n (k)\d j u n (]i.)}d 2 k, 
(44) 



where |u„(k)) = e~* k ' r |V'n(k)) is the periodic part of the 
Bloch wave function. In an insulator, the Hall conduc- 
tance, being quantized for any 2D plane in the 3D BZ, 
must be the same as we vary the out-of-plane component 
of k. Therefore, we can simply take this component of k 
to be zero. Using Eq.(14), we have 

e ij (d l u n (k)\d j u n (k)) = e ij (5 j es-u n (k)|5 i e)su n (k)jj45) 
= ei i (9 ? -u„(-k)|9iu„(-k)) 
= -eij(diU n (-k)\djU n (-k). 

Substitute this equation into Eq.(44) taking the other 
k m ^i,j = 0, we have 



4 = 0. 



(46) 



Eq.(46) states that, for any combination of i,j, there is 
no obstruction to defining a smooth gauge for all occu- 
pied states in the MBZ, making Eq.(43) a meaningful 
definition. 

However, this does not mean that So defined in Eq.(43) 
is a topological invariant, because we have not proved 
that it is gauge invariant. To see this, we perform a 
gauge transform that is continuous in MBZ: 



K(k)) 



m^occ 



U nm (k)\u m (k)) , 



(47) 



where U(k) is an arbitrary unitary matrix. The sewing 
matrix in the new basis becomes 



B mn \k) 






[/;„ m ,(-k)S m ,„,(k)[/,l,„(k). (48) 



From Eq.(48) we have 

det[B'(k)] = dct[[/*(-k)] det[t/ f (k)] det[B(k)] (49) 

= (dct[[/*(k)]) 2 det[£(k)], 
Pi[B'(k mv )} = dct[[/*(k m .„)]Pfp(k m ,)]. 

If in defining the complex square root func- 
tion y/z, we choose the branch cut such that 
if the argument of z is within [0, 27r), we have 

that dct [U* (k mv )}/ ^Jdet 2 [U* (k mv )] ±1 if 

arg[det[U*(k inv )}] e [0,ir) and arg[det[U*(k inv )}} £ 
[tt, 2it), respectively. If for odd number of ki nv we have 
arg[dct[U* (ki nv )\\ £ [tt, 2it), the ^-number defined in 
Eq.(43) changes its value and is therefore gauge variant. 
In a ID 0s-symmctric insulator, the only B-TRIM 
is k — 0. If we can find a gauge choice with 
arg[det[U*(0)]} £ [tt, 2tt), the Z2 number becomes gauge 
variant. Consider a constant gauge transform: U mn = 
^mn(l — 2<5 m i). By choosing this gauge we simply rede- 
fine the wave functions of the first band by multiplying a 



factor of — 1 while the wave functions of the other bands 
remain the same. Therefore, we have det[E/*] = — 1. Un- 
der this new gauge the Z 2 -number changes its value and 
is a gauge variant quantity, and hence no topological ©s- 
symmetric insulators exist in ID. 

If one considers a 2D 0s-symmetric insulator, there 
are three possible lattice translation vectors (D's): Di = 
af/2,D 2 = af/2,D 3 = (af + af )/2, where af 2 are 
the basis vectors of the magnetic superlattice. The ques- 
tion comes down to if we can find a gauge satisfying 
ar-g[det[E/*(k im) )]] € [7r,27r) at one of the two B-TRIM. 
When D = Di, there are two B-TRIM located at Ki = 
and K2 = bf /2. Consider the gauge transform 

C/„ m (k) = 6 mn (l + (exp(ik • af ) - l)S ml ). (50) 

In this gauge, the wave function of the first band is multi- 
plied by a factor of cxp(ik • af) while the wave functions 
of the other bands remain unchanged 34 . Using Eq.(49) 
we easily obtain det[£'(Ki)] = 1, det[B'(K 2 )] = -1. 
This gauge choice changes the value of the ^-number. 
For the other two possible D vectors, one can choose 
the gauge C/ mrl (k) = <5 mn (l + (exp(ik • af ) - l)8 m i) and 
U mn (k) = 8 mn {\ + (cxp(ik • (af - af )) - l)8 ml ) and 
reach the same conclusion. Therefore the Z 2 -number as 
defined in Eq.(43) is gauge variant also in 2D. 

The gauge variance in ID and 2D can be physically 
understood by applying the established theory classifying 
TRS insulators. To see this, one notices that the equation 
k • D = defines ad— 1-dimensional subsystem of the 
d-dimcnsional AFM insulator, in which we have 



B; 



[Q s ,H d _ 1 ] = 0, 
exp(2k-D) = -1. 



(51) 



These two equations are the defining properties of the 
symplectic class 29 (class All). Hence, the subsystem 
equals a spinful TRS insulator in d — 1-dimension. Since 
the topological classification of 0D and ID symplectic 
insulators is trivial, the ID and 2D Og-symmetric insu- 
lators must also be trivial. Therefore the Z 2 -number as 
defined in Eq.(43) must be gauge variant in ID and 2D, 
and no ID and 2D 0s-symmctric TFs exist. 

In the same light, the gauge invariance of the Z 2 - 
numbcr in 3D 0s-symmctric AFM insulators is easy to 
understand. The Hamiltonian constrained to the 2D 
plane satisfying k// • D = is exactly the same as a 2D 
TRS insulator with time-reversal symmetry replaced by 
0S, the combined symmetry, for which the ^-invariant 
given by Fu and Kane can be written down 



(-1)' 



Pfp(kl) 



n 

k // Jdet[B(k(i v )} 



(52) 



There are eight TRIM and four B-TRIM in any 3D insu- 
lator, and due to exp(2ik^ • D) = 1, all four k^'s are 
B-TRIM in the 3D MBZ. As a result, the Z 2 quantity 
defined in Eq.(52) is exactly the same as the one defined 
in Eq.(43) for 3D AFM insulators. Since in 2D TRS 



insulators the Z 2 quantity defined in Eq.(52) is gauge- 
invariant, we conclude that the Z 2 quantity defined in 
Eq.(43) is also gauge invariant. 

While we have proved that all ©s-symmetric AFM 
insulators have Kramers' degeneracy at B-TRIM, only 
in 3D can we define the Z 2 topological invariant. How- 
ever, when there are additional symmetries beyond 65 
in the system, new invariants may be defined in lower 
dimensions 12,13,15 . For example, when spatial inversion 
symmetry is present in the system, we can have well- 
defined Z 2 invariants in all dimensions. To see this, as- 
sume P to be the inversion symmetry operator in the 
single-particle Hilbert space. The commutation relation 
between the inversion operator and an arbitrary transla- 
tion operator is 



T r P = PT- r 



(53) 



and the commutation relation between the inversion op- 
erator and time-reversal operator is 



OP = PQ. 



(54) 



Therefore we have the commutation relation between 65 
and P as 



PSs = PQT D 

= ef_ D p 
= e s PT 2 D. 



(55) 



At B-TRIM, we have T 2D |^ n (k ml ,)) = \i>n(k inv )) . This 
means in the subspacc spanned by the two degenerate 
states at any B-TRIM, the inversion operator and 0s 
commute. If one state has parity (inversion eigenvalue) 
+1 (—1), the other must have the same parity. This 
allows us to define the following Z 2 invariant £o : 



(-l)fr = 



n 



Qn\' i inv)j (,5b) 



ki„„ GB-TRIM, nSocc/2 

where ( n is the parity of the n-th occupied band and 
the multiplication over occ/2 means only one state in a 
Kramers' pair is chosen. This definition applies to 1,2 
and 3D and can be extended using the same definition 
to any dimensions. These insulators, being stabilized by 
the added inversion symmetry, will, however, not exhibit 
gapless boundary modes in the energy spectrum but only 
in the entanglement spectrum. 



V. GAPLESS BOUNDARY MODES 

In 2D and 3D TRS insulators, there are odd num- 
ber of Dirac boundary modes in Z 2 non-trivial insula- 
tors and even number (including zero) of Dirac boundary 
modes in Z 2 trivial insulators, i.e., on the surface Bril- 
louin zone (SBZ) of a 3D TRS TI, there must be at least 
one Dirac node at one of the TRIM. At that Dirac point 
a Kramers' pair protected by time-reversal symmetry is 
located. Breaking TRS, in general, can gap the Dirac 



10 



point. For example, when the Fermi energy is close to 
the Dirac point energy, a magnetic field in a 2D TRS 
TI localizes the edge modes 35 ; and in a 3D TRS TI, ap- 
plying an FM Zeeman field on the surface opens gap at 
each Dirac point and makes the surface a 2D quantum 
anomalous Hall insulator, with each gapped Dirac node 
contributing ±e 2 //i 36 . 

When an AFM field is applied to a TRS TI, however, 
the Dirac point at a B-TRIM is still doubly degenerate 
if the system, with the open surfaces, still possesses ©s- 
symmetry. This requirement not only implies that the 
bulk preserves 0s-symmctry, but also implies that the 
surface termination does not violate that symmetry. This 
essentially puts constraints on how the system should be 
terminated in real space. 

We want to find out what types of surface terminations 
can preserve the Dirac nodes, i.e., the 0s-symmetry. 
First we remember that although we have been using D 
as if it were a unique vector, in fact D is unique only up 
to a magnetic lattice translation, i.e., any D' = D + L A/ 
is also a translation vector that inverts all spins. We 
so far have made no distinction between them because 
in a periodic system, if s = OTd is a symmetry then 
Q' s = 0Td' is also a symmetry and vice versa. For an 
open system, in order to preserve ©s, we need at least one 
D that (i) inverts all spins and (ii) is parallel to the sur- 
face. Rewrite Eq.(24) after expanding L M = ^V n^ 



This tight-binding model can be considered as rcprc- 



,i\/ 



D 



E («> + ?) 



M 



(57) 



where ni is an arbitrary integer and Xi is either zero or 
unity with at least one Xi being unity. (Here we assume 
that d is a dimension in which there is non-trivial topo- 
logical class for TRS insulators 10 .) If we use N to denote 
the normal direction of an open boundary, the protec- 
tion of surface states require the existence of a set of 
{ni, ...,nd} in Eq.(57) such that 



N-D = 0. 



(58) 



The above result may be understood from a more phys- 
ical point of view. When the condition Eq.(58) is met, 
all spins on the boundary change sign under an in-plane 
translation of D, and hence the total magnetization on 
the boundary is zero. To the lowest order perturba- 
tion, the total magnetization, i.e., ferromagnetism, on 
the boundary opens a Zeeman gap at a TRIM. No gap is 
opened to first order when the total magnetization van- 
ishes. On the other hand, if N is not perpendicular to 
any D, the total magnetization on the surface becomes 
finite and a gap is opened. The gap can be estimated 
by calculating the total magnetization within the decay 
length of the surface states. To be more concrete, let us 
consider a 2D AFM insulator, the TRS part of which is 
given by 



'Ho(k x ,k y ) = (m - cos^) - cos(k y ))T 
+ Y, sin(fci)(r< + 6r«). 



(59) 



(a) 



(b) 
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FIG. 4. (a) Different terminations of a 2D AFM insulator with 
propagation vector (tt,it), for which the coordinates of sites 
on the boundary satisfy y = nx. (b) The energy separation 
at the Dirac point of the system described by Eq.(59) with 
k z = 0, m = 1.5 and M = 0.2, as a function of the inverse 
slope of the cut for odd n. 

senting the 2D TRS TI of HgTe for momenta near T as 
1 < \m\ < 2 25 . The configuration of the AFM field is 
given by (see Fig. 4(a)) 



3(x,y) = Mzcos(irx)cos('Ky). 



(60) 



It is easy to see that this AFM insulator has 0s symme- 
try, with D = (af +af )/2 + L M , where a^f 2 = a 2 ± &1 . 
We choose several boundaries (edges) with the form: 



N n = -nai + a 2 . 



(61) 



i—x,y 



For n £ even, Eq.(58) is satisfied, by finding D = 
n/2(af - a^) + (af + a^)/2 = a x + na 2 ; but if 
n <G odd, Eq.(58) is unsatisfied and a gap is opened at 
the Dirac point. There are two ways to see the latter 
point: first, one can try to solve Eq.(58) substituting 
N„ with xi = x% = 1 and find that no integer solu- 
tion for rij's exists, or one simply notices that along the 
cut, the spins are all fcrromagnctically aligned, so it is 
impossible to find any translation along the edge that in- 
verts all spins. Along the edge, the distance between 
two nearest sites is \/n 2 + la, so the number of sites 
per unit length is p = l/\fn 2 + la ~ l/na and the dis- 
tance between the first and the second layer of atoms is 
d„ = a/\Jl + n 2 ~ a/n. The decay length of the edge 
mode at TRIM can be estimated as I ~ vp/^buik, where 
v p is the Fermi velocity at the Dirac point without AFM 
and Atuik is the bulk insulating gap. The total magne- 
tization close to the surface, or the ferromagnetic gap, is 
therefore 

AE n ex p(\ - cxp(-d/l) + cxp(-2d/0 - ...) (62) 
= p/(l + cxp(-d/Z)) 
oc 1/n, 

when n is large. By choosing the parameters m = 1.5, 
b = 0.01 and M = 0.2, we calculate the energy separation 
of the Dirac point at k = f in the edge Brillouin zone 
(EBZ). In Fig. 4(b), the energy separation is plotted ver- 
sus 1/n, and the linearity of the curve confirms Eq.(62). 
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This corroborates our intuition that when n is large, a 
surface with n £ even and n € odd cannot be physically 
distinguished. 

We have so far discussed whether the boundary Dirac 
modes are stable against a pcrturbativc AFM field. Does 
this hold beyond the regime of perturbation? Topologi- 
cal protection of gapless boundary modes in TI's means 
that: (1) the Dirac point is not split and (2) the boundary 
bands always lie in the bulk gap connecting the conduc- 
tion and the valence bands through spectral flow. The 
second point is important because otherwise one can adi- 
abatically push whole boundary band into the conduc- 
tion or the valence bands, making the system undistin- 
guishable from a trivial insulator. We will from now on 
focus on whether the boundary gapless modes can re- 
main within the bulk gap, or if they can be adiabatically 
pushed into the conduction or the valence bands. 

We begin by considering a 2D TRS TI subject to a 
perturbative AFM field that is Os-symmetric, in order 
to show from another perspective that AFM TI does not 
exist in 2D. We terminate the system and make an edge, 
the normal direction of which satisfies Eq.(58) for some 
D. Along the edge, a translation by D flips all spins, 
therefore we know the magnetic unit cell on the edge is 
of length 2D and the edge Brillouin zone (EBZ) extends 
from =£ to 2i}. Within the EBZ there are two TRIM: 
A"i = and K 2 = ^ . By the definition of A-TRIM and 
B-TRIM (see below Eq.(25)), we know that only K x is 
a B-TRIM. Therefore, without additional symmetries, in 
the EBZ, the edge bands must be doubly degenerate at 
A'i (zone center) and non-degenerate at K 2 (zone corner) . 
The non-degeneracy at zone corner is the key difference 
of this edge band dispersion from the dispersion of edge 
bands of all 2D TRS TI's and prevents the existence of 
spectral flow from the conduction band to the valence 
band. 

In Fig. 5(a), a schematic edge dispersion is plotted. 
There are two Dirac cones because the dispersions of both 
edges (left and right) are shown. The topology of a band 
structure remains unchanged as one arbitrarily distorts 
its shape without breaking any protected degeneracy. In 
our case, without a degeneracy at the zone corner, one 
can distort the dispersion from Fig. 5(a) to Fig. 5(b) and 
finally to Fig. 5(c) where a full gap is obtained, without 
closing the bulk gap. Since we have always preserved the 
degeneracies at the zone center, the fully gapped system 
in Fig. 5(c) is topologically undistinguishablc from the one 
in Fig. 5(a) where we have gapless edge modes. From this 
simple picture, we conclude that in a 2D AFM insulator, 
there is no protected gapless edge mode, although in the 
perturbation regime, the Dirac points of a TRS TI do 
not split and remain gapless. 

We explicitly realize the scenario sketched in Fig.5(a- 
c) in the 2D model given in Eq.(59), but with a magnetic 
structure 

3(x, y) = Mzcos(ny). (63) 

The system is cut along y-dircction and according to 



Eq.(58), the degeneracy at T in the EBZ is preserved. If 
we choose M = 0.2, the dispersion, plotted in Fig. 5(d), is 
similar to that of the edge mode of a TRS 2D TI. Then 
we increase the strength of the AFM field to M = 0.4 
(see Fig. 5(e)) and add different chemical potentials on 
the left and the right edges [II = — [1r = 0.5, the disper- 
sion becomes the one shown in Fig. 5(f), exhibiting a full 
gap between the conduction and the valence bands. 




FIG. 5. (a-c) Edge state dispersion under different param- 
eters without breaking the symmetries. It shows an adia- 
batic process, from (a) to (c), in which a full gap is open 
between the edge bands, (d-f) Edge and bulk dispersion of 
a 2D TRS TI described by Eq.(59) subject to AFM given 
by Eq.(63), in which the edge is along j/-axis. From (d) to 
(f), the bias potential is increased to separate the edge bands 
without either breaking the symmetry or closing the band 
gap. In (d), M = 0.2 (perturbative) and in (e,f), M = 0.4 
(non-pert urbative) . 

In 3D, consider a system with a surface having normal 
vector N. We can define the magnetic supercell, spanned 
by a^ 23 such that a^ 2 are parallel to the surface, i.e., 
perpendicular to N. (Generally speaking, such a choice 
breaks point group symmetries of the original lattice, but 
these symmetries are irrelevant in this paper.) Defining 
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bf 2 3 from aj u 2 3 , wc sec that the SBZ is spanned by bj " 2 , 
with three TRIM being f = 0, X = bf /2, Y" = bf/2 
andM = (bf +bf)/2. 

Since D is parallel to the surface, it must be equal, 
up to an in-plane supcrlattice vector, to one of the fol- 
lowing translation vectors: Di = a^/2, D2 = sl^ 1 1^ or 
D3 = (af 1 + a 2 v/ )/2. According to the definitions of A- 
TRIM and B-TRIM we find that the B-TRIM are f , Y 
when D = Di. Therefore, if we look at the band struc- 
ture along FY, there will be a protected surface mode 
always inside the bulk gap. On the other hand, when we 
look at the band structure along f X, we recover the situ- 
ation discussed for 2D AFM insulators, in which the sur- 
face band is doubly degenerate at T but non-degenerate 
at X, and can therefore be fully gapped by deformation. 
However, since along TY the surface mode always con- 
nects the conduction and the valence bands, the surface 
band as a whole is still topologically protected, and can 
never be pushed into the conduction or the valence band. 
The B-TRIM are f , X when D = D 2 , and f , M when 
D = D3. Similar statements on the surface modes can 
be made for these two cases by knowing the B-TRIM. 
To summarize the discussion of 3D AFM TI, we have 
proved that for an open surface that does not break the 
0s-symmetry, there is an odd number of surface bands 
crossing the band gap. The dispersion of these surface 
bands differ from those of 3D TRS TI in that along the 
direction having exp(2ik • D) = 1, the dispersion runs 
from the valence to the conduction band just like the 
surface band of a TRS TI, but along the other direction, 
the dispersion does not touch the conduction or the va- 
lence band. 

Below we use an explicit model to demonstrate the 
above statements. The TRS part of the 3D model is 
given by 



band; and the maximum energy decreases as n increases. 
A more experimentally accessible quantity than the en- 



Ho(k) = (m — y^ cos fcj)r 

i 

+ j^smkijTj +bTn 



(64) 



with parameters m = 2.5 and b — 0.01. This tight- 
binding model is equivalent to the one for Bi2Sc/Te3 in 
the continuum limit 28 , and on each open surface has a 
Dirac cone centered at T of the SBZ. The AFM field is 
given by 

2-kz 



3{z) = Mzcos ■ 



(65) 



We choose the yz-plane as the open surface. Based on 
Eq.(58), wc know that the Kramers' pair at T of the SBZ 
is preserved if n £ even. The dispersion on the SBZ 
for n = 2,4,6,8 are plotted in Fig.6(a-d) for M = 0.2. 
From these figures, we can see that the band dispersion 
along y-axis and that along z-axis are very different be- 
cause Y is a B-TRIM and Z an A-TRIM. Along y-axis, 
the band dispersion always cross from the valence band 
to the conduction band, while along z-axis, the disper- 
sion reaches a maximum energy below the conduction 





FIG. 6. The dispersion of the central surface band in the MBZ 
on surface in a 3D AFM TI with magnetic structure given by 
Eq.(65) with n = 2,4,6,8, respectively. The spectral weight 
for the first two layers of the surface is overplotted with the 
dispersion. 

ergy dispersion is the electron spectral weight that can 
be directly measured in an angle resolved photoemission 
spectroscopy 18 ' 37 ' 38 (ARPES), 

A(w,k) = ^2 Im / ^( c >i,ar(k)c^ Qr (k,i))cxp(iu;i), 

a,r,n 

(66) 

where n denotes the layer the electronic operator is acting 
on (because the translational symmetry has been bro- 
ken in the stacking direction). In reality, ARPES is a 
surface probe and therefore cannot penetrate an arbi- 
trary number of layers of atoms, and here we assume 
that it only probes the spectral weight of the surface 
rather than the bulk states. For a rough estimate, we 
only add up the contribution from the first two layers in 
Eq. (66). In Fig. 6, the surface spectral weight is over- 
plotted with the band dispersion, where lighter/darker 
means higher/lower spectral weight. From these figures, 
we can see that for AFM TI with large wavevector (small 
n), e.g., n = 2, the spectral weight seen in experiment 
differs only minutely from what can be seen on the sur- 
face of a 3D TI, because the zone corner where the band 
bends backwards is actually filled with states distributed 
mainly in the bulk. As n-increases, the surface states 
begin to show some downward bending along the z-axis. 

In Fig. 7, we have plotted the spectral weight at certain 
energies, assuming that ARPES can only pick up spectral 
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FIG. 7. The k-resolved spectral weight on the open surface 
of a 3D AFM TI with a magnetic structure given in Eq.(65). 
The spectral weight is plotted for E — 0.3, 0.5, 0.7 in Row 1,2 
and 3 respectively. The first column corresponds to M — 0, 
that is, a TRS TI; the second one corresponds to M = 0.2 
and n — 4; the third column corresponds to M = 0.2 and 
n — 8. The color bar shows the distribution of a state on the 
first two layers of the surface where unity means completely 
bound to the surface. 



weight from the first two layers. From these figures we 
can see that if there is no AFM, the equal energy contour 
is always closed, but when M ^ 0, at very low energy it 
is still closed but becomes open in the z-direction as one 
increases the energy. The critical energy at which the FS 
becomes open decreases as n increases (or the wavevector 
decreases). As opposed to the equal energy contours in 
the dispersion plotted in Fig. 6 which are always closed, 
the intensity plot of the spectral weight can show open 
arcs. This is because the states far away from the Dirac 
point are mostly bulk states having small spectral weight 
on the surface. 

It should be noted that the statement on the existence 
of protected boundary modes runs completely in parallel 
to the statement concerning whether a Z 2 -invariant can 
be defined. In a 2D AFM insulator, one cannot define 
a i^-invariant because the quantity defined in Eq.(43) is 
gauge variant in 2D, while in this section we have proved 
that it does not have protected edge modes that exhibit 
spectral flow, unlike 2D TRS TI, either. In a 3D AFM 
insulator, one can define a ^-invariant by Eq.(43) and 
at the same time, when it is .Z2-non-trivial, there is an 
odd number of protected surface modes inside the bulk 
gap. 



there has been no report of any TI that is AFM or- 
dered. On the other hand, in principle, we may use 
AFM substrate to induce AFM ordering in a TI thin 
film. For example, substrates of under-doped cupratcs 
and the parent compounds of iron-based superconduc- 
tors have strong AFM orderings and can induce AFM 
in the TI thin film with propagation vector (-7T, 7r, 0) and 
(-7T, 0,0), respectively 32,39 . For more complex AFM pat- 
terns one can use substrates with spiral antiferromag- 
netism or sinusoidal antiferromagnetism. In practice, 
we need the lattice matching between the substrate and 
the TI. For example, AFM substrates with tetragonal 
(cuprates) and orthogonal (iron-based superconductors) 
lattices in general do not match the rhombohedral lattice 
of Bi2Se/Te3. But based on the fact that the lattice con- 
stants in the 122-family of iron-based superconductors 
(~ 4A) are close to the lattice constant of the recently 
suggested TI of SmBg 40 " 42 (~ 4.13A), we speculate that 
one may use the AFM state of 122 iron-based supercon- 
ductors to induce stripe-like AFM in a thin film of SmBg. 
If the lattices do not match, we need them to be at least 
commensurate with each other. In that case, at the in- 
terface, one needs to define a new superlatticc, such that 
under translations by the new superlatticc vectors, both 
lattices are invariant. Then we can calculate the propaga- 
tion vectors in the new basis and check if ©s-symmetry is 
preserved. If it is, from the theory we know that the AFM 
preserves the Dirac degeneracy of the surface states, but 
makes the dispersion anisotropic as shown in Fig. 7. 

Part of the discussion for AFM insulators can be easily 
extended to superconductors with coexisting AFM order- 
ings. Specially, suppose we consider the superconductors 
with TRS and spin-orbital coupling (class Dili) , and with 
an intrinsic or induced AFM coexisting with SC. If the 
system is ©s-symmetric, the d — 1-dimcnsional subsys- 
tem defined by k • D =0 belongs to class Dili, where the 
TRS is replaced by ©s which for this subsystem satis- 
fies 0| = — 1. Therefore a d-dimensional AFM SC with 
@s symmetry has the same topological classification as 
a d — 1-dimensional TRS superconductors (class Dili). 
From the general classification in Rcf.9, 2D and ID Dili 
superconductors have Z% classification; therefore, 3D and 
2D ©s-symmetric superconductors also have Z^ classifi- 
cation. The state with coexisting AFM and SC may be 
realized by inducing an AFM in a 3D TSC thin film; 
for 3D TSC materials, copper doped Bi 2 Se/Te 3 43 and 
indium doped SnTe 44 have been proposed. We have as- 
sumed that no symmetry exists other than @s, but in 
general AFM materials can have non-trivial magnetic 
groups, which possibly lead to much richer topological 
classifications beyond the scope of the current work. 



VI. DISCUSSION 



What are the possible realizations of an AFM TIs? 
In Rcf.[22], GdBiPi was suggested as a TI with intrin- 
sic AFM ordering due to the presence possibly strong 
and frustrating spin exchange coupling. Experimentally, 



VII. CONCLUSION 

In this paper, we carry out a general study of the role 
of commensurate antiferromagnetism in topological in- 
sulators, generalizing the previous results of Ref.22. The 
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Hamiltonian of the electrons is modeled by a general TRS 
tight-binding model subject to a space dependent peri- 
odic Zeeman field having zero average (zero ferromag- 
netism). It is found that when there exists a lattice trans- 
lation D that reverses the AFM field at every site, the 
Kramers' pairs at half of the TRIM are unbroken by the 
AFM field, while at the other half of the TRIM the eigen- 
values arc gencrically non-dcgcncratc. The existence of 
D is guaranteed when the Fourier modes of the AFM 
field satisfies the following condition: for any set of m in- 
tegers Zi,..., m , if J2i=l,..., m z iQi = G > tnen Z)t=l,...,m*« 

must be even, where Q,'s are the propagation vectors of 
the AFM field. If this condition is not met, the energy 
splitting of the Kramers' pair can be estimated by know- 
ing the propagation wavevectors of the AFM field. Fur- 
thermore, we show that the Kramers' pairs at half of the 
TRIM make possible the definition of a new ^-quantity 
similar to the one defined for TRS insulators. However, 
the new ^-index is only meaningful in 3D, but not in 
ID or 2D (as it becomes a gauge variant quantity). The 
surface manifestation of the new Z2-index is the presence 
of gapless Dirac modes in 3D AFM TIs. There are three 
differences between the surface modes of a AFM TI and 
those of a TRS TI: (i) to see the surface states, in AFM 
TI the surface of termination must preserve the combined 
symmetry 0g while in TRS TI the system can be termi- 
nated at any surface; (ii) in a AFM TI, the dispersion of 
the surface band only connects the conduction and va- 



lence bands in one direction, but is separated from the 
bulk bands along the other direction, while in a TRS TI 
the surface band connects the conduction and the valence 
bands along any direction; (iii) the spectral weight of the 
surface states at a certain in-gap energy makes a closed 
loop in a TRS TI, but appears as disconnected arcs in an 
AFM TI at energies away from the Dirac point. We dis- 
cuss the materials in which the theory can be applied and 
also how the theory can be easily extended to the case 
of topological superconductors with coexisting AFM or- 
ders. Since in our work we have not assumed any other 
symmetry or any specific magnetic structure, the theory 
applies to general insulators with commensurate AFM. 
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Appendix A: The equivalence between the two 

sufficient and necessary conditions of the 

Os-symmetry 

In the text, we derived two conditions for the AFM 
to have Os-symmctry from different physical perspec- 
tives. First, 0s-symmetry is equivalent to the existence 
of a lattice vector D, the translation by which inverts all 
spins. From considering AFM as perturbation, however, 
0s-symmetry means that for any given set of m integers 
zi,2,...,m, if J2i z iQi = G) then J2i z % = even. These two 
conditions are asserted to be mathematically equivalent 
to each other; but the equivalence is not made transpar- 
ent in the text and is proved in this Appendix. This 
Appendix is purely mathematical and may be skipped 
unless one needs to be convinced of the assertion by a 
formal proof. 

First we formulate the statement in a more mathe- 



matical fashion. Suppose the commensurate AFM has 
m propagation vectors, each having d components (in a 
rf-dimcnsional space). We can put these vectors into a 
d x m rational matrix, where qtj — a^ • Qj/27r is the i-th 
coefficient of Qj expanded in terms of b^. The lattice 
vector D can be expanded in a^'s, D = Xa=i d^* a «- 
In these symbols, the two conditions can be put into the 
following two statements: 
For a given matrix q, 

S A : 3D = {D u ...,D d )eZ d ,s.t. 2D T q€ (2Z + l) m ; 
Sb'- V z G Z m s.t. qz e Z d , there is J2i=i m z i~ even - 
To understand why the satisfying Sa is equivalent to the 
existence of D inverting all spins, simply notice that if 
Sa is true, there is D • Qj = odd * ir for any j. 

We will follow the steps, detailed in three subsections, 
to prove that the statements Sa and Sb are equivalent: 
1. Define a set of linear transformations on any given 
q, q' = RqC, where R is an integer d x d matrix and 
C a m x m matrix with integer entries. Then we prove 
that if Sa,b(<1) = 1, then Sa,b(q') — 1 and vice versa. 
(Here by Sa,b(q) = 0/1, we mean that Sa,b is false/true 
for matrix q.) In other words, Sa,b are invariant under 
these transformations. 2. Prove that using the above 
transformations, any given q can be transformed into one 
of the basic forms, to be defined later. 3. Prove that for 
any matrix go of a basic form, Sa,b{<1o) = 1 or Sa,b(qo) = 
0. There are five basic forms: B\ = Odxm', Ri satisfies 
(-82)1.7 = Pij/^ij where p\j ^ and fc lj = fcy for any 
j,j'; B^s has at least one column that is zero; B[ is a 
diagonal matrix with first i c elements being in the form 
qa = Pu/2 kii where pa 5^ 0, and all the other elements 
are zero; B' 2 is a block-diagonal matrix with the upper 
block a diagonal matrix of dimension if — 1 satisfying 
qa = Pii/^ kii where pa ^ and the lower block in the 
form of B^ ■ 



1. Linear transformations that leave Sa,b invariant 

Below we define seven transformations on q that leave 
SA,B{q) invariant: 
(i) interchange two rows in q; 
(ii) interchange two columns in q; 
(iii) subtract an integer riy from g^-; 
(iv) multiply the i ^ 1-th row by an integer and add it 
to the first row; 

(v) multiply the first row by an odd integer; 
(vi) multiply the first column by an integer Ci satisfying 
12i=i m Ci = even and c\ = 0, and add all resulted 
columns to the first column; and 
(vii) multiply the first column by an odd integer. 

By using (i,ii), we can easily extend the definitions of 
(iv,v,vi,vii) to 

(iv') multiply the i 7^ j-th row by an integer and add it 
to the j-th row; 

(v ! ) multiply any row by an odd integer; 
(vi') multiply the i-th column by an integer Ci satisfying 
12i=i m c i = even and Cj = 0, and add all resulted 
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the reduction, after which every q^ can be expressed as 
Pij/sij with Sij > pij. If all elements are zero after the 
step, we have reached a basic form B\ = Odxm- 

Step two: If this is not the case, then there is at least 
one qij jt= remaining. (From this point we no longer 
differentiate q and q' , as we have shown that Sa,b is 
invariant under the transformations listed.) Note that 
every s^ can be written as a product of 2 kij and some odd 
integers, where k^ > 0. Using transformations-^', vii') 
we can eliminate all the odd integers from Sij . After this 
we again use transformation- (hi) to eliminate all integer 
entries, (qi/s whose Sij only has odd integers are put to 
zero.) Perform the procedure for each g^ whose Sij has 
an odd factor. This is the second step of reduction. The 
reduction stops if after this step q = B\ . 

Step three: If after the second step q is still nonzero, 
we have Sij = 2 ij for every nonzero element of qij . Find 
the largest kij and send the corresponding element q^ 
to <7n using transformations- (i,ii). In the next step, we 
first use transformation-(iii,iv') to eliminate all the other 
elements in the first column. This can be done as follows: 
(1) for any qn = pn/2 kil ^ 0, multiply the first row by 
2*;ii -feii anc } a ^ y. ^ ^g j.^jj row . ^2) now we have 

Qii = (9ii + < ?ii)/2 fcil , so one can cancel at least a factor 

of 2 in the numerator and denominator of qn , after which 

where CV, = Sij for i ^= 1 and C\\ — r G odd. First we the new s' n = 2 kil where k' n < kn and (3) repeat (1,2) 



columns to the j-th column; and 

(vii') multiply the j-th column by an odd integer. 

Except for the first two transformations, which mean a 
permutation of the spatial dimensions and a permutation 
of the propagation vectors, the physical meaning of the 
other transformations is not obvious. But the physical 
meaning is not our concern in the proof, and we only 
need that under any of these transformations, Sa,b{q) 
is invariant for any q. The proofs of the invariance are 
straightforward for all transformations, yet we prove the 
invariance under (vi) and (vii) explicitly here. 

The transformation in (vi) is expressed as q' = qC, 
where CV, = 5ij + riSn satisfies J] r« = even and r\ = 0. 
First we prove that Sa{q) = 1 — > SaW) = 1- From 
(2D T q)i G odd, we have Y. J {2D T q) 3 C jl = (2D T q)i for 

i^ 1 and E j (2D T q) j C jl = (2D T q) 1 +J2 i n(2D T q) l G 
odd, i.e., (2D T q')i G odd. Then we prove that Ssiq) = 
1 -> Ss(q') = 1 by exclusion. If Ssiq') = 0, there 
is z' such that q'z 1 = qCz' G Z d and Yl z 'i "= °dd. 
Defining z = Gz', we have qz G Z d and J2(Cz')i = 
J2 z'i + Y^riz'i G odd, so Sb{q) = 0. Third, we notice 
that C _1 can be obtained by changing r t — > —ri, so C _1 
also belongs to this class of transformations. Therefore, 
Sa,bW) = 1 — > Sa,b(q) = 1 follows automatically. 
The transformation in (vii) is expressed as q' = qC, 



prove S A {q) = 1 -> S A (q') = 1. Since (2D T q) l G odd, 
and {2D T q') l = {2D T qC)i = (2D T q) l for i ^ 1 and 
(2£> T g')i - {2D T qC) 1 = r(2i? T < Z ) 1 , we have (2D T q') l G 
odd. Second we prove £U(<z') = 1 — >• Sa(q) = 1. Since 
(2D T q')i G odd, we have [2(rZ?) T g']i G r*odd; then there 
is [2(rL») T g] l = [2(r J D) T g'C'- 1 ] 4 e r * odd * CT. 1 G odd. 
Third, we prove Sb(si) = 1 — > Sb(q') = 1 by exclusion. If 
S B (g') = 0, then there is z' such that 2gV = 2gCz' G Z d 

and J2 z i ^ °dd, so X)(^ z ')i G °^> arL( i z = ^V is an 
integer vector that satisfies 2qz G Z d and J2 Zi G odd. 
Finally we prove Sb(<i') = 1 — ^ Sb(q) = 1 by exclusion. 
If Ssio) = 0, then there is z such that 2qz G Z d and 
£) 2j € odd. From 2qz G Z d we have 2qC(C~ 1 rz) G 
Z d and ^(C^Vz)., G odd. So we have z' = C~ x rz G 
Z m such that 2q'z' G Z d and £ z'j G odd. We have 
shown that under linear transformation described in (vii) , 
Sa,b{q) is invariant. 

As a final remark, the invariant transformations in- 
troduced here do not form a group, because the inverse 
of transformations- (v,vii) involves dividing a row or col- 
umn by some integer, which do not leave .S^b invariant 
in general. 



2. Reduction to basic forms 

Here we show that using the above transformations, 
every g-matrix can be reduced to one of the 'basic forms'. 

Step one: First we stress that for commensurate AFM, 
every element q^ must be a rational number. If any 
\Qij\ ^ 1: we can u se transformation- (hi) to send it to 
a value between —1 and +1. This is the first step of 



until the qn becomes an integer, then use transformation- 
(iii) to eliminate it. Next we try to eliminate qij^i in the 
first row. If there is any qij^i = 0, then we can use 
this clement together with qa to eliminate all the other 
elements in the first row through transformation-(vi'). 
Here the elimination process is almost the same as the 
one eliminating qu^o (after replacing rows by columns) 
with one catch: when k\j = feu, 2 kll ~ klJ = 1 G odd so 
subtracting the first column from the j-th column is not 
one of the invariant transformations, but since there is 
some qij 1 — 0, we can subtract the first row and the j\- 
th row from the j-th row. If q±j 7^ for any j, but there 
is some k\j 2 < k\\, we can use the first row to eliminate 
qij 2 , then use the first and the J2-th row to eliminate all 
the other elements in the first row. This is the third step 
in the reduction process. The elimination of the first row 
cannot proceed if and only if q\j ^ and kij — fen, that 
is, when all elements in the first row are nonzero and 
have the same denominator. A g-matrix in this form is 
defined to be of a basic type called B2 ■ 

Step four: Suppose after the previous step q ^ B2, 
then we have q\j = qn ~ for i, j ^ 1 and qn ^ 0. The 
g-matrix is now block-diagonalized into two blocks, the 
upper left block having only one element g" = gn, and 
the lower right block containing the rest (g')y' = Qi+ij+i- 
For q l , we repeat all the above three steps, with three 
possible outcomes: qi = B\, qi G B 2 or q l is block- 
diagonalized with its first row and column eliminated 
except q n . We define B[ as the set of matrices hav- 
ing qn 7^ for i < i c > where < i c < d and all 
other elements zero; and define B' 2 as the set of block- 
diagonalized matrices whose upper blocks are diagonal 
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matrices of dimension if — 1 > and lower blocks are of 
basic form B 2 . We see that if we repeat the reduction 
process finite times, we have the following three possibil- 
ities: q is reduced to basic form B[; q is reduced to basic 
form B 2 ; when d < m, qu 7^ for i < d and all other 
elements are zero, this is of basic form B3 because there 
is at least one zero column. 

We have shown that under a series of the invariant 
transformations introduced, any g-matrix is reduced to 
one of the five basic forms: B 1: B 2 , B[, B' 2 and B 3 . 



3. Proving Sa = Sb for all basic forms 

In the part, we finish the proof by showing that for any 
g-matrix in a basic form, denoted by go, SU(go) = Sb(<Zo)- 

If qo = B\, it is trivial to show that Sa(qo) = Sb(qo) = 
0, i.e., both statements are false. Physically, this basic 
form corresponds to an FM state, which of course does 
not exhibit 0s-symmetry. 

If go £ B[, we show that Sa(qo) = Ssiqo) = 0, i.e., 
both statements are false. To see this, simply notice 
that in go there is at least one zero column, namely the 
jo-th column, therefore (2D T go)j = ^ odd for any 
D £ Z d and £U(<Zo) = 0. On the other hand, define a 
vector zq £ Z m such that (zq)j = Sjj , then we have 
q z = £ Z d and £(,2o)j = 1 £ odd, so S B (qo) = 0. 
Similarly, we can show that for go £ -B3, we also have 
SU(go) = Ssiqo) = 0. Physically, this basic form cor- 
responds to ferrimagnetism, which has modulating local 
magnetization with a nonzero average. 

If go £ B 2 , we show that SU(go) = Sb^o) = 1- We 
choose D £ Z d : D t = 2 fc "" 1 for i < if and A = for 
i>if, and we have 

2D T g = (pn,P22, -,Pifif,Pifi f +i, -,Pi f m,)- (Al) 

Since all p^'s appearing in Eq.(Al) are odd, SU(<Zo) = 1- 
For any z £ Z m , if q z £ Z d , there must be z, = 2 fc " 
for i < if, and ^7=1 m z jPifi *= even -, from which 
we know V" ■ ■ „, z; must be even, hence Y^ z, £ even. 

L — /l=l f ,...,771 ^ L — i b 

The proof for Sa = Sb is very similar if g £ B 2 , because 
B 2 can be seen as a special case of B' 2 with if = m. 



We have proved that for go in any basic form, Sa (<Zo) = 
Sb(<1o)- Combining the results from the previous two 
subsections, Sa = Sb for any g-matrix. 



Before closing the Appendix, we gives a few exam- 
ples how we reduce a commensurate AFM, i.e., a ma- 
trix q, to its corresponding basic form go. We as- 
sume a 2D system in these examples. For the first ex- 
ample, consider an AFM with one propagation vector 
Q = (bi +b 2 )/3, so q = (1/3, 1/3) T . By transformation- 
(vii), we have g — > (1, 1) T and by transformation- (iii), 
we have go = (0,0) T £ B\. In the second exam- 
ple, consider an AFM with two propagation vectors 
Qi = bi/24 + b 2 /4 and Q 2 = 3bi/24 + 5b 2 /4, i.e., 

A/24 3/24\ 
q = I 1 ,. 1 1. I . First use transformation- (v) to ob- 
tain g — > \Jia J I a ) 1 then multiply the first row by 

—2 and add it to the second row (transformation-(iv')), 
A/8 3/8^ 



which becomes go = [ L , /„ ) £ B 2 . In the third 

example, consider an AFM with three propagation vec- 
tors Qi = bi/4, Q 2 = b 2 /2 and Q 3 = bi/2, i.e., q = 
1/4 1/2^ 
1/2 
(vi'), multiplying the first column by —2 and add it to 

1/4 0\ 
1/2 Oj E Bz - 
In the next example, Qi = bi/4 and Q 2 = b 2 /3, i.e., 
A/4 
' ^ 1/3; 
(vii') followed by transformation- (iii) becomes go 
"1/4 0^ 
y 

/ 1/2 N 
with Qi = b 2 /6 and Q 2 = bi/2, i.e., g = L , ' 

First we use transformation- (ii) to interchanged the two 
columns, and then use transformation- (v') to obtain 

_ A/2 0\ 
qo - { 1/2 ) G S 2' 



To simplify, just use transformation- 
;he first colum 
the third column, resulting in go 
nple, Qi = bi/< 

, which after applying transformation- 

>y transformation-(iii) becomes go = 

£ B' \. For the last example, consider an AFM 
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